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Abstract 

In this paper, I study final state QCD radiative corrections to off- 
shell single top production via decaying W at hadron collider energies. 
The tree level Breit-Wigner distribution of the produced top invariant 
mass is distorted by final state QCD radiation, while the peak posi- 
tion remains unchanged. The exact one loop QCD calculation and the 
narrow width approximation agree in predicting the cross section, the 
hadronic transverse energy distribution and the bottom-lepton invari- 
ant mass distribution. 
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1 Introduction 



The discovery of the top quark at CDF and DO |T||, has opened a new era 
of measurements in top-quark physics. Now, the properties of the top quark 
can be directly investigated, not only inferred from their effects in radiative 
corrections. 

At hadron colliders, the dominant production mechanisms are, of course, the 
tt channels 

qq — > tt 

gg^ti, (l) 

but single top quarks events are also present, such as 

q'g(W + g) -> qtb 

q' b — > qt 

q'q _> w* -> tb 

gb -> W~t . (2) 

The first two mechanisms are known as W-gluon processes 0, the third one 
as W* production || and the fourth one as Wt production Q. The cross 
sections in eq. (Q) are ordered according to their magnitude in pp collisions 
at y/s = 2 TeV for m t = 176 GeV §. 

Even with less expected events, single top production processes are important 
because they provide a consistency check on the measured parameters of the 
top quark in tt production. 

Radiative corrections to the processes in eq. (|J) and eq. (|J) are well 
known in the literature 0, but usually, performed in the narrow width ap- 
proximation, in which production and decay of the top are treated indepen- 
dently. A check on the validity of this approach is still missing. Since, in the 
narrow width approximation, diagrams connecting decaying products with 
the production process are missing, one especially expects deviations due to a 
non exact treatment of the gluonic radiation, which is an important quantity 
for the reconstruction of m t in tt events 0. Therefore a precise study of it, 
even in a simpler case, can give hints on its relevance in the main production 
mechanisms of eq. ([!]). 

For those reasons, I decided to perform a complete calculation of the fi- 
nal state QCD radiative corrections to W* single top production, taking into 
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account all the subsequent decays. Among all the others, the W* mechanism 
is interesting because possible new physics may introduce a high mass state 
(say particle V) to couple strongly with the tb system such that the produc- 
tion rate from q' q — > V — > tb can deviate from the standard model W* 
rate || ||. Therefore accurate predictions of the standard properties of the 
top in this channel are also important by themselves. 

The background QCD contribution is known ||, so I shall study here the 
QCD one loop corrections to the signal diagram of fig. 1. For this process, 
thanks to the color structure, initial and final state QCD corrections do 
not interfere and are separately gauge invariant, so that, in order not to 
obscure the effects I want to study, I decided, in this paper, to concentrate 
my attention on final state gluonic corrections. Initial state corrections are 
however very simple and a study including them will be treated elsewhere 
[0]. I chose the semi leptonic final state V\ l + bb because it is easier to detect 
experimentally. The extension to q' q bb is anyway trivial. In fact, diagrams 
with gluons connecting b or b with q' or q are killed, at the one loop level, by 
the color factor, so one is left with simple gluonic corrections for the W q' q 
vertex. 
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Figure 1: Tree level diagram for al u — > v\ l + b b via single top production. Here, 
and in the following figures, dashed lines denote W 's. 

2 The calculation 

The tree level diagram for the process is drawn in fig. 1, while in fig. 2 and 3 
I show the one loop virtual diagrams and the real bremsstrahlung. 
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Figure 2: Final state one loop QCD virtual diagrams. 




Figure 3: Real gluon emission. 

I computed the virtual corrections using standard Passarino-Veltman tech- 
niques |TIJ, with the help of the Symbolic Manipulation program FORM 



nj. I used dimensional regularization for ultraviolet, collinear and soft di- 
vergences. Furthermore, I kept everywhere complex masses for top and W, 
but I systematically neglected the bottom mass. 

An analytic approximation in n dimensions for the soft-collinear part of the 



real emission was built following ref. [12] and the cancellation of all diver- 
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gences performed analytically. 

Both the virtual and the real contributions have been computed applying 
helicity amplitudes methods JHJ, and the final expressions implemented in 
a Monte Carlo program [0], that uses the self-optimization techniques of 



ref. fll4| . Since those techniques are applied here, for the first time, in loop 
calculations, it may be useful to briefly discuss the adopted strategy. More 
details will be found in ref. 0. 

The problem here is the matching between hard and soft phase space. 
Schematically, the final result for the cross section o t ot (with any kind of 
cuts) can be written as a sum of four contributions 

&tot = o"o + cry + cr s (5) + <7 H {S) (3) 

where o"o is the lowest order result, ay the virtual contribution and as, an 
the soft and hard real radiation. 

The sum, ay + as($) does not contain soft and collinear singularities. On 
the other hand as (5) + a H (5) is independent on 5, where S is the separa- 



tion between soft and hard gluons (following again ref. |12], 5 in a cut on 
the invariant mass of g + b and g + b). The last statement is true only 
if an exact computation of as (5) is performed. Instead, what one usually 
does is computing as{o) for small 5. In such a limit, because of factorization 
properties (|l|, fll)]), very simple expressions are obtainable in terms of the 
born result multiplied by universal coefficients containing \og(5) and log 2 (ci). 
At this point, by numerically going to the limit 8 — > 0, one gets unbiased 
results. Of course, if S is too small, large numerical cancellations take place 
between as(S) and ajj{o), resulting in large errors in the Monte Carlo inte- 
gration. A good value for S can be usually found by numerically checking 
the independence on 5 of the results. 

For fixed 5 one would like to know how many Monte Carlo points have to be 
spent to separately compute all four contributions in eq. @ , mainly because 
usually the most time consuming part is ay, that contains loop diagrams. 
This is a typical problem that can be solved using the Multichannel self- 



optimizing approach of ref. Jl4|]. One starts with the same amount of points 
for all channels and, during the run, the Monte Carlo self-adjusts itself, so 
that afterwards one usually obtains a smaller percentage of the computational 
time spent in the evaluation of ay, which means a better estimate of atot in 
a shorter time. 
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channel 


percentage before opt. 


percentage after opt. 


1 


0.2 


0.0996 


2 


0.2 


0.5436 


3 


0.2 


0.1265 


4 


0.2 


0.2083 


5 


0.2 


0.0220 



Table 1: Percentage of the Monte Carlo points used for each channel in the 
evaluation of a to t before and after the self-optimization. The first three chan- 
nels take care of the peaking structure of an given by the Feynman diagrams 
in fig. 3, channel 4 refers to o~ + a$ an d channel 5 to ay . 

In table 1 I show a typical result of the self-optimization procedure. 

I made several checks on the final result. First of all, I verified that the 
CP invariance of the tree level current 

T a „ = u(b) j a (l - 7 5 )($+ +p l Vl + ^ fe + m t )7 M (l -j 5 )v(b) (4) 

remains after QCD loop corrections. Then, by numerically rescaling Tyy, T t 
and the cross section by the same amount, I checked the agreement between 
the Monte Carlo estimate of the total 0(a s ) tb on-shell cross section (amc) 
and (for example) the analytic result (a an) of ref. |16| (see table 2). I also 
tested, for small 5, the independence on S of the results. All numbers in this 
paper are obtained with 5 = 0.2 GeV 2 . 

A last comment is in order. Taking into account the widths of the de- 
caying particles gives rise to conceptual problems with respect to the gauge 
invariance. The correct gauge invariant prescription would be to compute 
the widths as the imaginary part of the one loop renormalized propagators 
and all set of loop diagrams necessary to restore gauge invariance. Since, 
in the process at hand, W and t decay via electroweak interactions, this 
would imply to include terms of the one loop 0(a) calculation at the tree 
level and part of the two loop 0(aa s ) corrections in the 0(a s ) contribution. 
Since I am interested here in 0(a s ) corrections, the neglected 0(a) terms 
are expected to be small, so I followed the naive prescription of considering 
everywhere constant complex masses. For the sake of consistency, when com- 
puting 7m(m|), I used the lowest order top width value T t = 1.6429 GeV for 
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(GeV) 


°mc (pb) 


o~an 


300 


0.09086 ± 0.00002 


0.09088 




0.10499 ± 0.00035 


0.10547 


600 


0.03508 ± 0.00001 


0.03507 




0.03784 ± 0.00027 


0.03774 


900 


0.01653 ± 0.00001 


0.01653 




0.01763 ± 0.00019 


0.01748 


1200 


0.00948 ± 0.00001 


0.00948 




0.00999 ± 0.00014 


0.00996 


1500 


0.00612 ± 0.00001 


0.00612 




0.00643 ± 0.00010 


0.00640 


1800 


0.00427 ± 0.00001 


0.00427 




0.00454 ± 0.00007 


0.00446 



Table 2: Comparison between the Monte Carlo total cross section, in the limit 
of vanishing widths, and the analytic on-shell calculation. No convolution 
with the parton densities has been performed. The first entry is the tree 
level result. In the second entry all final state QCD corrections are included. 
m t = 176 GeV and a s = 0.12. 

quantities at the tree level and the QCD corrected value Y t = 1.5117 GeV 
when including QCD corrections. 

3 Results 

In this section, I present some results for the process d u — > v\ l + bb obtained 
with the Monte Carlo of ref. ||. I chose to plot three useful distributions 
for measuring the top mass in pp collisions at yfs = 2 TeV, namely the total 
hadronic transverse energy (H T ), the invariant mass \J (pi+ + pb) 2 (my) and 

the "top mass distribution" J(p Vl +Pi+ +Pb) 2 (jnbiv)- Of course, due to the 
presence of an undetected neutrino, the last quantity is not going to be easy 
to reconstruct experimentally. However, m&z„ is of theoretical interest and 
directly measurable in the channel q' q bb. 
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I used the following input parameters and cuts 

a = 1/128, sin^ = 0.2224, a s = 0.1 
M w = 80.41 GeV, M t = 176 GeV 
T w = 2.1185 GeV 

Eriyi), E T (l + ), E T (b), E T (b) > 15 GeV 

\rj(l+)\, \ V (b)\, \rj(b)\ < 2, AR(bb)>0.7 } (5) 

together with the cone jet-definition algorithm of ref. [|17[] (with jet cone 
size R= 0.7) and the CTEQ2M parton densities of ref. IJI8"] . Two jets with 
b content are required to be present in the visible region defined by the 
above cuts, and no extra (gluonic) jets. In order to compare the full QCD 
calculation with the narrow width approximation, I also produced histograms 
for Ht and mu , in which W and top are put on-shell by numerically rescaling 
IV,i and the cross section by the same amount. 

In fig. 4 and 5, I show Ht and mu in the off-shell case and in the on-shell 
limit, including all final state gluonic corrections. For comparisons, I also 
show the tree level result. As one can see, the on-shell and off-shell distribu- 
tions are almost indistinguishable, therefore the narrow width approximation 
works very well. 

In fig 6. the off-shell top invariant mass distribution is shown with and 
without final state QCD corrections. QCD radiation is responsible for a 
distortion in the Breit-Wigner distribution: more events are produced in the 
left tail and less in the right side, but the position of the peak is essentially 
unchanged []. Comparing with the narrow width approximation is difficult 
for m b i u . One would be forced to impose a Breit-Wigner by hand. I did not 
try that. Fig. 6 already gives a quantitative prediction for the distortion 
of the one loop distribution with respect to the tree level result. A precise 
quantitative knowledge of this effect may be useful in fitting the experimental 
distributions, and estimating the systematic errors. 

In on-shell production, the cross section (with the cuts and input parameters 
of eq. (^|), including final state QCD corrections) is 0.02677 ± 0.00025 pb, 
while one gets 0.02681 ± 0.00024 pb in the off-shell case. That means that 
the narrow width approach can be safely used also in computing the cross 

1 By looking at fig. 4 and 5 one can recognize a similar distortion in Ht and mu as 
well. 
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section. The question of the total number of produced events is important 
when looking at the single top production rate in this channel for new physics 
searches. 

4 Conclusions 

I have performed a complete one loop calculation of the final state QCD 
corrections to the single top production process d u — > v\ l + b b in the off- 
shell case. Final state gluonic radiation is responsible for a distortion in the 
distributions useful for top mass reconstruction. The validity of the narrow 
width approximation is confirmed at the level of accuracy one naively expects, 
namely 0(T t /m t ). However, one should observe that, no one loop QCD di- 
agrams connecting initial and final states can contribute to the process at 
hand. On the other hand, such diagrams are present in the main produc- 
tion mechanisms of eq. ([]]). Therefore, in order to check the narrow width 
approach in that more general full off-shell QCD loop calculation for 

the it channels is still needed. 
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Figure 4: Total hadronic transverse energy for on- shell (dashed histogram) and 
off-shell (solid histogram) single top production. Final state QCD corrections 
included. The dotted line is the off-shell tree level result (a s = 0). 
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Figure 5: The histograms are the invariant mass distribution of l + + b for 
on-shell (dashed) and off-shell (solid) single top production, including final state 
gluonic corrections. The dotted line is the off-shell tree level result (a s = 0). 
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Figure 6: Invariant mass of l + + b + v\ for off-shell single top production, at the 
tree level (dotted histogram) and including final state QCD corrections (solid 

histogram). 
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